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Abstract. In this contribution we review a series of simple one dimensional lattice models that 
with an appropriate choice of parameters can account for various anomalous features of the be¬ 
haviour of complex systems such as water. In particular, we will focus on the presence of p — T 
fluid-solid coexistence lines with negative slope (i.e. solids that melt upon compression), solid 
phases less dense than the liquid phase, and the existence of temperatures of maximum density. 
We will see how a simple two-parameter model can reproduce the phase behaviour of a range 
of systems well known for their anomalous behaviour regarding the temperature and pressure 
dependence of properties such as density, diffusivity or viscosity. 


1. Introduction 

The singular properties of water and its essential role in our lives have encouraged an enormous 
amount of theoretical and experimental work since the first half of the twentieth century. But it has 
been the search for a possible second critical point [I] [2J [SJ and the careful study of the many facets of 
the anomalous behaviour of water what has been a key point for research in the last decades [H [1,6, 7j. 
By anomalous behaviour in water we mean well known features, such as the negative slope of the 
p-T liquid-solid coexistence line, which in common words means that the solid phase (ice) will melt 
upon compression (see the phase diagram of Figure[l]). Together with this, water is known to present 
a temperature of maximum density (TMD), which at atmospheric pressure is approximately 4°C, 
this implies that there exists a region of anomalous behaviour in which water expands when cooled 
down. At the same time one finds a region in which diffusivity increases when pressure is increased 
(dynamic anomaly). When plotted in a T — p diagram the various anomalous regions organize in 
a cascade of anomalies[8]. But not only water exhibits this type of anomalies, other substances 
as diverse as, PP HU], Bi, C, Si, Ge, Ga, Te,p, silica [T2], or germanium oxide|13j have pressure- 
temperature coexistence curves with negative slope, present liquid-liquid equilibrium, or have been 
shown by simulation to exhibit dynamic anomalies. All have in common the presence of stable solid 
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Figure 1. Experimental phase diagram of water. Roman numerals designate dif¬ 
ferent ice phases. Data taken from Refs. muni nsj. 


phases with relatively low coordination numbers (3-5, lower that those of the corresponding liquid 
phases). 

Interestingly, it has become clear that such anomalous behaviour can also be reproduced by 
model systems interacting via spherically symmetric potentials)!?). In particular, the Hemmer-Stell 
potential)!!], devised with the aim of describing simple systems that can exhibit multiple transitions, 
was shown by Jagla llfll to share with water a good number of anomalies. In this regard, Yan et. 
al[17), where the first to describe in this model a cascade of anomalies which is qualitatively similar 
to that of water (see Figure [2]). Together with the Hemmer-Stell’s (or ramp) potential, a variety 
of interaction models have shown to be able to display peculiar behaviour, such as the Lennard- 
Jones+Gaussian (LJG) [20]; or the square-well/square shoulder (SWSS) [2] among others. It has 
been argued that only the continuous version of the SWSS potential [2T] yields a density anomaly, 
whereas other features, such as the presence of a liquid-liquid equilibrium are reproduced by the 
SWSS potential. These models are all characterized by the presence of two repulsive ranges of 
interaction. There is however an additional simple potential model that shares a good number of 
anomalies with these systems, namely the Gaussian core model|22]. Since this potential, being finite 
at zero separation, allows for complete overlap of a pair of interacting particles, it is a class of its 
own. One may think that together with a first repulsive range defined by the width of the Gaussian, 
a second range of interaction would correspond to full overlap (or the maximum of the repulsive 
potential), thereby casting this model into the class of two repulsive range models. 

Additionally, some lattice gas models have also been investigated to search for thermodynamic 
and structural anomalies characterized by of both isotropic (24] 125] and orientational dependent 
interactions [26 [ 127] . In this contribution we will review the simplest models that can reproduce 
a water-like phase behaviour and at least some of the above mentioned anomalies. Thus we will 
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Figure 2. Cascade of anomalies for a two-short range repulsive potential (Hemmer- 
Stell) as studied by Yan et al.[T7|. Qe defines the orientational order parameter, t 
the translational order parameter (see [53]). DM designates the curve of diffusivity 
maxima and TMD the line of temperatures of density maxima. 


focus on one-dimensional lattice models. In this regard, the simplest system that can give a fluid- 
solid transition (order-disorder) in one dimension is the antiferromagnetic Ising model with two 
sublattices, in which the antiferromagnetic state (ordered) is mapped into a solid in lattice gas 
language. The model, studied in detail by Hpye [28] yields a second order transition near the critical 
point. We will review the features of this model in the next section. Then in Section III we will 
see how it can be modified so as to reproduce a first order fluid-solid transition with a solid phase 
less dense than its fluid counterpart, following previous work by the authors[55]. In Section IV we 
will see how the addition of a second repulsive range brings the phase diagram closer to that of a 
particular class of anomalous systems such as phosphorous [30]. The two-parameter model can be 
thereafter tuned to yield a water-like phase diagram |31l . A purely repulsive version of the model will 
also be shown to reproduce the presence of density maximum, bringing much closer the analogy to 
water behaviour. The careful study of these and related simple yet rich physical models has been a 
distinctive characteristic of Prof. Johan Hpye’s scientific career, representing in many cases relevant 
landmarks in Statistical Physics. 

2. First step: antiferromagnetic models with two and three sublattices. 

Back in the early seventies Hpye [38] studied closely the order-disorder phase transitions of an 
antiferromagnetic Ising model with two sublattices, in which the antiferromagnetic state (ordered) 
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Figure 3. Qualitative phase behaviour of the one-dimensional Ising model with 
nearest-neighbour repulsion and an infinitely long ranged staggered mean-field act¬ 
ing on even numbered neighbours [251. 


can be mapped into a solid phase. The model is endowed with competing short range and very long 
range interactions, being the latter of Kac type, and acting only on even numbered spins, i.e. 


<p(i-j) = { — a T e xp[—7|* — j\\, 
where a is the integrated long-range interaction, i.e. 


if i — j is even 
if * — j is odd 


a = 


1 OO 

2 51 ¥>(*)> 


( 2 . 1 ) 


( 2 . 2 ) 


and 7 —> 0 being and inverse range parameter. So this long range interaction does not necessarily 
encourage ferromagnetic ordering. The full Hamiltonian can be written as 

N N N 

H = SiSi +1 +^2<p{i- j)siSj - H 0 ^2 Si- (2.3) 

i— 1 i<j i 

with J > 0 being the antiferromagnetic nearest neighbour (NN) coupling constant, N is the number 
of spins or lattice sites, s, : = ±1 are the spin variables, H 0 is the external field, and we consider 
periodic boundary conditions, i.e. Sj = s^+i- This can be shown to be one of the simplest models 
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that can yield a fluid-solid transition. However, in Ref. [28] , H 0 ye, showed that such a system leads 
to a second order transition in a region close to the critical point, which for temperatures below a 
certain transition temperature, turns into the desired first order fluid-solid transition (see Figure 


i- 
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Figure 4. Phase diagram of a lattice gas with infinite NN repulsion and and three 
sublattices with a staggered mean field. 


In order to cure this deficiency, Hpye and Lomba [25], extended the model incorporating three 
sublattices, i.e. we considered an infinitely long ranged antiferromagnetic interaction that acts on 
each third spin, i.e. 




3 o 7 exp[— 7 |i — j'l], if i — j = 3n + 1, 3n + 2 

0, if i — j = 3 n 


(2.4) 


where n is integer, and again the limit 7 —» 0 is considered with a > 0. The model thus constructed 
would have however an unwanted solid-solid transition in the middle of the order-disorder transition, 
due to the symmetry of the Ising model when the critical point is located at zero magnetic field 
or zero magnetization. The presence of a nearest neighbour hard core interaction will break the 
symmetry with respect to the zero field, and so remove this unwanted transition. This model can be 
solved using a mean field treatment to handle the long range staggered interaction and the transfer 
matrix method to deal with the NN antiferromagnetic coupling(32] [28]. As shown in [29] one can 
now get a system that yields a fluid-solid transition, in which the solid phase melts again upon 
compression (see Figure [4]). In order to simplify the solution, the antiferromagnetic NN coupling 
was considered in the limit J —> 00 , by which then the external field Hq must be adjusted so that 
Hq + 2 J remains finite and spins can be turned up and down by thermal fluctuations. The model 
is equivalent to a lattice gas in which the particles occupy two sites. The phase diagram of Figure 
[4] indicates that we are now in a situation closer to a water-like behaviour: we have a model in 
which the solid phase melts upon compression and with a p — T coexistence curve with negative 
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Figure 5. T-p (left) and p-T (right) phase diagrams of a one-dimensional lat¬ 
tice model with NNN interactions and a staggered field. Dotted curves indicate 
metastable equilibria. 


slope within certain range of temperatures, i.e., pressure decreases when temperature raises. This 
represents a clearly “anomalous” behaviour. 

3. One step beyond: two sublattices and next-to-nearest neighbour interaction 

How can one improve the previous model ? On one hand we have to account for the vapour- 
liquid transition. This in principle can be achieved separating the mean field contribution to the free 
energy into a staggered term (that favours antiferromagnetic long range ordering) and an attractive 
term. This has to be complemented by the incorporation of a next-to-nearest neighbour interaction 
(NNN), which will enable to tune the location of fluid-solid equilibrium with respect to the gas-liquid 
critical point. 

The system so constructed was studied by the author in collaboration with Hpye in [80] . and 
solved using the procedure indicated above and which will be roughly described below. The model’s 
Hamiltonian can be written as 

NNN N 

H = SjSj+1 + K ^2 SiSi +2 + ^2 vi * ~ j ) s i s i ~ H 0 ^2 Si ■ (3.1) 

i — 1 i — 1 i<j i 

Again the J > 0 is NN coupling, and as in [29] we let J —> oo with Ho adjusted so that Hq + 2 J 
remains finite, and periodic boundary conditions are imposed, but we have now a NNN interaction 
with a coupling constant K > 0. 
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The lattice can be split into groups of cells made up of three sublattices in accordance with the 
staggered interaction (2.4). The effective fields, (acting on spins of sublattice i) can be expressed 
in terms of the mean field interaction (2.41 and the external magnetic field as 

Hi 
H 2 

H 3 = H 0 -a(mi+m 2 ), (3.2) 


= Hq — a(m 2 + m 3 ) 
= Hq — a(mi + m 3 ) 
= H 0 — a(mi + m 2 ), 


where m, is the magnetization per particle of sublattice i. Now, we define mi = m + 2u and 
m 2 = m 3 = m — u, where m and u are the effective and staggered magnetizations respectively. Thus 
the following effective field, H e , and effective staggered field H s , can also be defined by 


H e = l -(H 1 +H 2 + H 3 ) = HQ-2am 

H s = - H 2 - H 3 ); u = \{2mi - m 2 - m 3 ). 

3 6 


These satisfy 


Hq = H e + 2 am 
He = 2au 


(3.3) 

(3.4) 

(3.5) 

The equivalency to the lattice gas gives a number density p = (to + l)/2. The reduced temperature 
will be defined in terms of the mean field coupling as T* = 1 /(/3a), with /3 = 1/ksT, ks the 
Boltzmann constant, and T the absolute temperature as usual. The corresponding reduced field and 
NNN coupling constant will be h* = Hq/cl and K* = K/a. 

Now, in order to account for the vapour-liquid transition the effect of the uniform mean held term 
2 am can be tuned, while retaining the staggered mean Held contribution, 2au. One thus ends up 
with 

H e = Hq — 2 aam (3-6) 

where a will be a coupling term that can switch from uniform mean Held attraction (a < 0) to 
repulsion (a > 0). Thus a/1 represents the additional uniform mean field interaction mentioned 
below Eq. ( |2.4| . 

The partition function, A, is given by the eigenvalue equation of the transfer matrix (Eqs. (23)-(25) 
in Ref. [30]), namely 

A 3 + a 2 A 2 + a\A + oo = 0 (3.7) 

where A = A 3 , and 


a 2 

ai 

a 0 


D(2e~ v + e 2v )p\ 


-[ D 3 p “ 3 
-~[D 2 {e~ 2v + 2 e v ){p~ 2 - p 2 ) + p ~ 6 ] 
D 3 p- 3 (p~ 2 -p 2 ) 3 , 


(3.8) 


with 


D 


= exp (~/3H e ) 

v = 13 H Sb ■ p = exp(/3K). (3.9) 

The magnetization to and staggered magnetization u result from the differentiation of the eigen¬ 
value equation by the method used in Ref. [32]. provided A has been calculated solving Eq. (3.71. 
Thus, one gets 


to = — 


1 3 
1 y, 

3 J 

n—0 


' n A n /N 


(3.10) 






Figure 6. Phase diagram of phosphorous. Data taken from Ref.[9]. 


3 3 

- E < An / N ’ N = tl 


n—1 


n—1 


with 


dan 


n d(/3H e ) 

, _ da n 

~ d((3H Se ) ’' 

Explicit expressions can be found in Ref. |30j. 

With this, the free energy G per spin [251 ED] is given by 

— (3G = In A + /3aam 2 — /3au 2 . 


(3.11) 


(3.12) 


(3.13) 


and the corresponding lattice gas pressure p reads 

/3 p = —)3G + f3H + /3K + /3 aa. (3.14) 

We have thus all the ingredients to study the phase behaviour of our model. To that aim, we 
know that the phase equilibrium conditions determine that the field, h*, and the spin free energy, 
(3G, stay the same in both phases at equilibrium. Thus one has to solve, 

f3G(u = 0, m d ; T*) = /3G(u a , m 0 ; T*) 
h*(u = 0,m d -,T*) = h*(u a , m 0 \ T*) 


(3.15) 

(3.16) 
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Figure 7. T -p (left) and p-T (right) phase diagrams of a one-dimensional lattice 
model with NNN interactions and a staggered field with a choice of mean field 
attraction (a) and NNN repulsion(AT/a) to mimic the phase diagram of water. 


where the subscripts d and 0 denote the disordered and ordered phases respectively. Additionally, 
md and h*, and m D , u 0 and h* are connected via Eqs.(3.10) and (3.11). Details concerning the 
explicit numerical solution of these equations can be found in Refs. m and [30] . 

After a careful search over the parameter space of K/a and a, one finds that the equations (3.16) 
yield the phase diagram for K/a = 2.5 and a = —8 as shown in Figure[5] Here both the T — p and 
T — p phase diagrams are shown (p is evaluated from Eq. ( 3.14[ )). Clearly this phase diagram departs 
from that of water (Figure [l]) in the fact that the solid-liquid equilibrium curve starts out from the 
triple point with a positive instead of a negative slope. At higher densities the L-S curve reaches 
a liquidus point and then a region of anomalous behaviour (i.e. negative slope in the T — p curve) 
appears. As a whole, this diagram bears some resemblance to that of phosphorous [9] depicted in 
Figure [6] Also in Figure [5] one observes the presence of a metastable liquid-liquid equilibrium, a 
feature that has been argued to be found in undercooled water (though much debated) [I, 6]. On the 
other hand, in the case of phosphorous, a first order liquid-liquid transition has been characterized 
by means of diffraction experiments in the region of high pressures and temperatures [33l llOj . 


4. A WATER-LIKE MODEL 

Now, the question is whether the model parameters can be tuned so as to reproduce the behaviour 
of water. To that aim, it seems clear from Figure [5] that the liquidus point has to be shifted into the 
vapour-liquid equilibrium curve. The LV curve has to be made somewhat larger in that respect so 
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Figure 8. Temperature of maximum density curve of a one-dimensional lattice 
model with NNN interactions and a staggered field with purely repulsive interac¬ 
tions. The region of anomalous behaviour is indicated by a shaded area. 


that the VL critical temperature is well above the liquidus temperature. We found that a choice of 
a = —15 (a substantial increase in the dispersive interactions) and K/a = 2.75 (a slight increase in 
the NNN repulsion), leads to a phase diagram in qualitative agreement with that of water, as can 
be seen when comparing Figures [l] and [7] One can observe in Figure [7] that now the LS equilibrium 
curve starts at the triple point with negative slope, as expected for a water like model. The change in 
density at the triple point is, on the other hand, close to one per cent, far away from the experimental 
Ap/p ss 10% in real water. If one is to bring this to a better quantitative agreement with “real” 
water, both a and very specially K/a must be fine tuned. This repulsive parameter is a key quantity 
since being a NNN repulsion plays an essential role in the determination of the low density solid 
phases. Work on this and other refinements is currently in progress. 

We have seen then that a tuning of our two parameter model switches from phosphorous-like to 
water-like behaviour. In order to understand the key differences between both systems, one must 
analyze the nature of bonding in the different phases of water and phosphorous. In the case of water, 
both the solid and the liquid phases are dominated by the presence of relatively strong hydrogen 
bonds. On the other hand, low density liquid phosphorous is a molecular liquid of P 4 tetrahedra 
which interact via weak van der Waals forces, and the solid is formed by a layered low density 
orthorhombic phase[5], which upon isothermal compression melts back into a high density liquid 
with entangled chains and clusters of covalently bonded P atoms [33]. It is evident that attractive 
forces are much more important in liquid water that in the molecular P 4 liquid. This explains why 
in order to go from a P-like to H 2 O-like phase diagram one has to substantially increase the strength 
of the attractive interaction (which rule the vapour-liquid equilibrium) from a = —8 to a = —15. 
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As mentioned, changes in K/a should make possible the fine tuning of the density change at the 
triple point. 

Now, a question that remains to be analyzed is whether this simple model is able to reproduce the 
existence of density maxima for certain temperatures (at constant pressure), which is a characteristic 
feature of water and related systems. Using Eq. (3.14) to evaluate p along isochorous curves, from 
elementary thermodynamics we know that those points that fulfill 


(s^), = 0al8Of " lflI1 (^)„ = 0 ' 

i.e., these state points correspond to density maxima along isobars. If one performs such a cal¬ 
culation with the parameters found for our water-like model, the resulting TMD curve lies in the 
metastable/unstable region. We have thus analyzed whether a corresponding purely repulsive model 
lacking a LV equilibrium can capture this feature by performing the calculations with a = 0. The 
result of these calculations is plotted in Figure [8] As we see the curves present clear minima which 
display a TMD curve that is a boundary of a thermodynamically anomalous region, in which the 
fluid expands upon cooling. This is actually the feature that was being sought. It remains to be 
seen whether a sensible choice of parameters can make compatible the existence of a TMD curve in 
a thermodynamically stable region with a water-like phase diagram, thus bringing our simple model 
to a much closer agreement with the experimental behaviour. 

In summary, we have seen how one can construct a family of one dimensional lattice models that 
can be tuned to reproduce the behaviour of fairly complex materials such a phosphorous or water. 
This simple models illustrate key physical phenomena enabling a more clear understanding of the 
underlying physics of the anomalies found in the thermodynamics of not so simple materials. They 
also help expanding our knowledge of how their anomalous behaviour is originated which will in turn 
help in the development of more complex models that can reproduce more closely the experimental 
behaviour. 
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